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We discuss phenomenology in a new TeV scale model which would explain neutrino 
oscillation, dark matter, and baryon asymmetry of the Universe simultaneously by the 
dynamics of the extended Higgs sector and TeV-scale right-handed neutrinos. Tiny 
neutrino masses are generated at the three-loop level due to the exact Z2 symmetry, 
by which the stability of the dark matter candidate is guaranteed. The model provides 
various discriminative predictions in Higgs phenomenology, which can be tested at the 
Large Hadron Collider and the International Linear Collider. 

1 Introduction 

In spite of the success of the Standard Model (SM) for elementary particles, it is widely 
understood that a new model beyond the SM must be considered to explain the phenomena 
such as tiny neutrino masses and their mixing pQ, the nature of dark matter (DM) [5] and 
baryon asymmetry of the Universe [3] . 

We here discuss the model in which these problems would be simultaneously explained 
by the TeV-scale physics [3]. Tiny neutrino masses are generated at the three-loop level 
due to an exact discrete symmetry, by which tree-level Yukawa couplings of neutrinos are 
prohibited. The lightest neutral odd state under the discrete symmetry is a candidate of 
DM. Baryon asymmetry can also be generated at the electroweak phase transition (EWPT) 
by additional CP violating phases in the Higgs sector [5]. In this framework, a successful 
model can be made without contradiction of the current data. 

The original idea of generating tiny neutrino masses via the radiative effect has been 
proposed by Zee [5]. The extension with a TeV-scale right-handed (RH) neutrino has been 
discussed in Ref. [TJ, where neutrino masses are generated at the three- loop level due to 
the exact Z2 parity, and the Z2-odd RH neutrino is a candidate of DM. This has been 
extended with two RH neutrinos to describe the neutrino data |S]. Several models with 
adding baryogenesis have been considered in Ref. [9] . The following advantages would be in 
the present model [4] : (a) all mass scales are at most at the TeV scale without large hierarchy, 

(b) physics for generating neutrino masses is connected with that for DM and baryogenesis, 

(c) the model parameters are strongly constrained by the current data, so that the model 
provides discriminative predictions which can be tested at future experiments. 

In the following, we first explain the basic properties of the model, and discuss its phe- 
nomenology, in particular that at the International Linear Collider (ILC). 
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2 Model 



Two scalar isospin doublets with hypcrchargc 1/2 ($1 and $2), charged singlet fields (S ), 
a real scalar singlet (77) and two generation isospin-singlet RH neutrinos (N R with a = 1, 2) 
are introduced in our model [4] . We impose an exact Z 2 symmetry to generate tiny neutrino 
masses at the three- loop level, which we refer as Z 2 . We assign Z 2 -odd charge to S ,± , 77 and 
iVjj, while ordinary gauge fields, quarks and leptons and Higgs doublets are Z 2 even. In order 
to avoid the flavor changing neutral current in a natural way, we impose another (softly- 
broken) discrete symmetry (Z 2 ) [JU]. We employ so called Type-X Yukawa interaction [TT] . 
where Z 2 charges are assigned such that only $1 couples to leptons whereas $2 does to 

quarks [HI nam; 

C Y = ^ e T$ie R -y u $® 2 u R -y d fi% 2 d R + \i.c., (1) 

where Q l (L l ) is the ordinary i-th generation left-handed (LH) quark (lepton) doublet, and 
u l R and d l R (e R ) are RH-singlct up- and down- type quarks (charged leptons), respectively 
We summarize the particle properties under Z 2 and Z 2 in Table [TJ 

The Yukawa coupling in Eq. ([TJ is different from that in the minimal supersymmetric SM 
(MSSM) [T5J. In addition to the usual potential of the two Higgs doublet model (THDM) 
with the Z 2 parity and that of the i?2-odd scalars, we have the interaction terms between 
Z 2 -even and -odd scalars: 

2 / 2 \ 2 

A»t = -£(p.|*«lW+*«l*«l a y) ~ E {k e<*(* c a ?*bS-ri + h.c}, (2) 

o=l ^ ' a,b=l 

where e a b is the anti-symmetric tensor with e\ 2 = 1. The mass term and the interaction for 
N R are given by 




Although the CP violating phase in the Lagrangian is crucial for successful baryogenesis at 
the EWPT [5J, it does not much affect the following discussions. Thus, we neglect it for 
simplicity. We later give a comment on the case with the non-zero CP- violating phase. 

As Z 2 is exact, the even and odd fields cannot mix. Mass matrices for the Z 2 -even scalars 
are diagonalized as in the usual THDM by the mixing angles a and /3, where a diagonalizcs 
the CP-even states, and tan/3 = (& 2 )/(§i) [TS]. The Z 2 even physical states are two CP- 
even [h and H), a CP-odd (A) and charged (H^) states. We here define h and H such that 
h is always the SM-likc Higgs boson when sin(/3 — a) = 1. 
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Table 1: Particle properties under the discrete symmetries. 
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Figure 1: The diagrams for generating tiny neutrino masses. 
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0.0021 


7.8xl0~ i2 



Table 2: Values of hf for m H± (m s ±) = 100 (400) GeV, m v = 50 GcV, m N i R = m N 2 =3.0 
TcV for the normal hierarchy. For Set A (B), ntan/3 = 29 (34) and U e3 = (0.14). 
Predictions on the branching ratio of \i — >• ej are also shown. 

3 Neutrino Mass, Dark Matter, and Strongly lst-Order Phase Tran- 
sition 

The LH neutrino mass matrix is generated by the three-loop diagrams in Fig. [T] The 
absence of lower order loop contributions is guaranteed by Z2- and e l R play a crucial role 
to connect LH neutrinos with the one- loop sub-diagram by the Z 2 -odd states. We obtain 

2 

M H = ^ C % F ( m H±> m S±, m N^ m v)> ( 4 ) 
a=l 

where C% = 4k 2 tan 2 /3(yf i M /if )(yf™hf) with y™ = V2m e Jv and v ~ 246 GcV. The factor 
of the three-loop integral function F{m H ±,m s ±,m N includes the suppression factor 

of l/(16-7r 2 ) 3 , whose typical size is O(10 4 )cV. Magnitudes of Ktan/3 as well as F determine 
the universal scale of M y -, whereas variation of hf (i = e, r) reproduces the mixing 
pattern indicated by the neutrino data [1]. 

Under the natural requirement h" ~ O(l): an d taking the /i — > ej search results into 
account [TB], we hnd that m Na ~ 0(1) TeV, m H± < 0(100) GeV, /ctan/3 > O(10), and 
m s ± being several times 100 GeV. On the other hand, the LEP direct search results indicate 
m H± (and rn s ±) > 100 GeV [Tj. In addition, with the LEP precision measurement for the 
p parameter, possible values uniquely turn out to be m H± ~ m H (or m A ) ~ 100 GeV for 
sin(/3 — a) ~ 1. Thanks to the Yukawa coupling in Eq. (JJ), such a light is not excluded 
by the b — > s"f data |17j . Since we cannot avoid to include the hierarchy among yf M , we only 
require hfyi ~ 0(y e ) ~ 10~ 5 for values of hf. Our model turns out to prefer the normal 
hierarchy scenario. Several sets for hf are shown in Table [2] with the predictions on the 
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Figure 2: [Left figure] The relic abundance of rj. [Right figure] The region of strong first 
order EWPT. Deviations from the SM value in the hhh coupling are also shown. 



branching ratio of /i — > ej assuming the normal hierarchjQ . 

The lightest Z 2 -odd particle is stable and can be a candidate of DM if it is neutral. In our 
model, must be heavy, so that the DM candidate is identified as r\. When r\ is lighter than 
the W boson, r\ dominantly annihilates into bb and t + t~ via tree-level s-channel Higgs (h 
and H) exchange diagrams, and into 77 via one-loop diagrams. From their summed thermal 
averaged annihilation rate (av), the relic mass density fl n h 2 is evaluated. Fig.[2jLeft) shows 
il 71 h 2 as a function of m, ; . Strong annihilation can be seen near 50 GeV ~ m H /2 (60 GeV 
~ m/j/2) due to the resonance of H (h) mediation. The data (f^DM^ 2 — 0.11 [5]) indicate 
that m v is around 40-65 GeV. 

The model satisfies the necessary conditions for baryogenesis [3]. Especially, departure 
from thermal equilibrium can be realized by the strong first order EWPT. The free energy 
is given at a high temperature T by 

V eff [<p,T\ = D(T 2 - T V - ET^ + + (5) 

where <p is the order parameter. A large value of the coefficient E is crucial for the strong 
first order EWPT with keeping mh < 120 GeV. For sufficient sphalcron decoupling in the 
broken phase, it is required that [18 



<Pc(„2E_ 



> 1, (6) 



where tp c 0) and T c are the critical values of (p and T at the EWPT. In Fig. [5jRight), 
the allowed region under the condition of Eq. © is shown. The condition is satisfied when 
m s± > 350 GeV for m A > 100 GeV, m h ~ 120 GeV, m H ~ m H± (~ M) ~ 100 GeV and 
sin(/3 — a) ~ 1, where M represents the soft-breaking mass of extra Higgs bosons for Z2 [I]. 



a The predictions for fi — > e-y shown here are corrected ones from those in Ref. pp. 
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4 Phenomenology 



A successful scenario which can simultaneously solve the above three issues under the data [TJ 
□H [17] would be 



sin(/3 — a) ~ 1, 
toa > 0(100) GeV, 



Ktan/3 ~ 30, 



400 GeV, 



m h = 120 GeV, 



0(100) GeV, 
- 3 TeV. 



(7) 



This is realized without assuming unnatural 
hierarchy among the couplings. All the masses 
are between 0(100) GeV and 0(1) TeV. The 
discriminative properties of this scenario are in 
order: 

(I) h is the SM-like Higgs boson, but decays 
into 7777 when < nih/2. The branching ratio 
is about 30% for m v ~ 43 GeV and tan/3 = 10. 
This is related to the DM abundance, so that 
our DM scenario is testable at the CERN Large 
Hadron Collider (LHC) and the ILC by search- 
ing the missing decay of h. Furthermore, rj is po- 
tentially detectable by direct DM searches [19], 
because rj can scatter with nuclei via the scalar 
exchange [20] . 



root(s)=500GeV 
root(s)=800GeV 
root(s)=1000GeV 




Figure 3: The production cross section of 
e+e- -> HA. 



(II) For successful baryogenesis, the hhh cou- 
pling has to deviate from the SM value by more than 10-20 % 
be tested at the ILC l2"2l 



(see Fig. which can 



(III) H (or A) can predominantly decay into 
t + t~ instead of bb for tan/3 > 2 because of 
the Type-X Yukawa interaction. For exam- 
ple, we have B(H(A) -> r+r") ~ 100 % and 



0.3% form, = 



= 130 



.S+ 



Nr 



s- 



.s- 



Figure 4: Feynman diagrams for the pro- 

s+s-. 



B(H(A) -> n+v 

GeV, sin(/3 - a) = 1 and tan /3 = 10. The sce- 
nario with light H ± and H (or A) can be directly 
tested at the LHC via pp -> IV* -)• i/F ± and 
7li7 ± [2^ , and also — > HA. Their signals 
are four lepton states £~£ + t ± i' and £~£ + t + t~ , 
where £ represents /i and r [11]. At the ILC, the 

process e + e~ — > if A would be useful to discriminate the model from the other new physics 
candidates. In Fig. [3] the production rate of the e + e~ — > HA is shown for m A = m H . 
For y/s = 500 GeV, about 17,000 (110) of the t+t-t+t- (/i+/j-r + T-) events are then 
produced from the signal for = m# = 130 GeV QT], while about 60 (0) events are in 
the MSSM for the similar parameter set. The main back ground comes from ZZ production 
(about 400 fb), which is expected to be easily reduced by appropriate kinematic cuts. 

(IV) The physics of ^-odd charged singlet is important to distinguish this model from 
the other models. At the LHC, they arc produced in pair via the Drell-Yuan process [24] . 
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The cross section amounts to 0.5 fb for 
m s± = 400 GeV at y/s = 14 TcV, so that more 
than a hundred of the S + S~ events are produced 
for the integrated luminosity 300 fb -1 . The pro- 
duced bosons decay as 5' ± — > H r], and _ff ± 



mainly decay into r 



The signal would be a 
with a large missing 
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Figure 
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high-energy hadron pair 
transverse momentum. 

The charged singlet scalar bosons 5 ± in our 
model can also be better studied at the ILC via 
e+e" S + S~ shown in Fig. H The total cross 
sections are shown as a function of m s ± for y/s in 
Fig. [5] The other relevant parameters are taken 
as m. Tl = m» r2 =3 TeV and hi = hi = 2.0. 

"r "r 

Both the contributions from the s-channel gauge 
boson (7 and Z) mediation and the t-channel 
RH neutrino mediation are included in the cal- 
culation. The total cross section can amount to 

about 200 fb for m s± = 400 GeV at yfs = 1 TcV due to the contributions of the t-channel 
RH neutrino-mediation diagrams with 0(1) coupling constants . The signal would be a 
number of energetic tau lepton pairs with large missing energies. Although several processes 
such as e + e~ — > W + W~ and e + e~ — > H + H~ can give backgrounds for this final state, we 
expect that the signal events can be separated by kinematic cuts. 

Finally, there is a further advantage in testing our model at the e~e~ collision option of 
the ILC, where the dimension five operators l~ t~ S + S + , which appear in the sub-diagram 
of the three-loop induced masses of neutrinos in our model, can be directly measured. The 
production cross section for e~e~ — > S~S~ [t-channel mediation: sec Fig. [6] is given by 



300 
m s (GeV) 



Production cross sections 
— > S + S~ via the s-channel 
gauge boson (7 and Z) mediation (dotted 
curve), the t-channel RH-neutrino (iVJ|) 
mediation (solid curve), and both contri- 
butions (dashed curve) for y/s = 300, 500 
and 1000 GeV. 



cr(e 



dt- 



1 



1287TS 



a = l 



1 



t — m 



1 



(8) 



Due to the structure of our model that the tiny neu- 
trino masses are generated at the three-loop level, 
the magnitudes of h% (a = 1,2) are of 0(1), by 
which the cross section becomes very large. Fur- 
thermore, thanks to the Majorana nature of the t- 
cahnnel diagram we obtain much larger cross sec- 



sion when m 2 Na 3> s. Fig. [7] shows the production 
cross sections for e~e~ — > S~S~ via the t-channel 
RH-ncutrino. The cross section can be as large as 



30 



pb for to 

-m Nl 



s± = 400 GeV for ^s e - e - = 1 TeV, 
3 TeV and hi=h 2 = 2.0. The back- 



grounds are expected to be much less than the e H 
collision. 



s- 



:cn r 



he**"- 
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Figure 6: Feynman diagram for the 
processes of e~e~ — > S~S~ . 
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We emphasize that a combined study for these processes would be an important test 
for our model, in which neutrino masses are generated at the three-loop level by the 
symmetry and the TeV-scale RH neutrino^. 

In the other radiative seesaw models in which the neutrino masses are induced at the 
one-loop level with RH neutrinos, the corresponding coupling constants to our couplings 
are necessarily one or two orders of magnitude smaller to satisfy the neutrino data, so that 
the cross section of the t-channcl RH neutrino mediation processes are small due to the 
suppression factor (h" ) 4 . 

(V) The couplings hf cause lepton flavor viola- 
tion such as fj. — > e"f which would provide infor- 
mation on rriN^ at future experiments. 

Finally, we comment on the case with the 
CP violating phases. Our model includes the 
THDM, so that the same discussion can be ap- 
plied in evaluation of baryon number at the 
EWPT [5|. The mass spectrum would be 
changed to some extent, but most of the fea- 
tures discussed above should be conserved with 
a little modification. 

5 Summary 

We have discussed the model with the extended 
Higgs sector and TeV-scale RH neutrinos, which 
would explain neutrino mass and mixing, DM 
and baryon asymmetry by the TeV scale physics. It gives specific predictions on the collider 
phenomenology. In particular, the predictions on the Higgs physics are completely different 
from those in the MSSM, so that the model can be distinguished at the LHC and also at 
the ILC. 

References 

[1] W. M. Yao, et al., [Particle Data Group] J. Phys. G 33 (2006) 1. 

[2] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 180 (2009) 330. 

[3] A. D. Sakharov, Pisma Zh. Eksp. Tcor. Fiz. 5, 32 (1967). 

[4] M. Aoki, S. Kanemura and O. Scto, Phys. Rev. Lett. 102, 051805 (2009); Phys. Rev. D 80, 033007 
(2009). 

[5] J. M. Clinc, K. Kainulaincn and A. P. Vischer, Phys. Rev. D 54, 2451 (1996); L. Fromme, S. J. Huber 
and M. Seniuch, JHEP 0611, 038 (2006). 

[6] A. Zee, Phys. Lett. B 93, 389 (1980) [Erratum-ibid. B 95, 461 (1980)]; A. Zee, Phys. Lett. B 161, 141 
(1985). 

[7] L. M. Krauss, S. Nasri and M. Trodden, Phys. Rev. D 67, 085002 (2003). 
[8] K. Cheung and O. Seto, Phys. Rev. D 69, 113009 (2004). 

b Unlike our model, in the model in Ref. [7J, the coupling constants corresponding to our h™ arc small 
and instead those to are O(l), so that its Majorana structure is not easy to test at e~e~ collisions. 



















Root(s)=1000GeV^\ 






500 \ 




300 


h,=h, =2.0 






m N H '= m N H 2 =3 T sV 



100 200 300 400 500 

m s (GeV) 

Figure 7: Production cross sections for 
e~e~ — > S~S~ via the t-channel RH- 
neutrino (iVg) mediation for s/s = 300, 
500 and 1000 GcV. 



7 



[9] E. Ma, Phys. Rev. D 73, 077301 (2006); J. Kubo, E. Ma and D. Suematsu, Phys. Lett. B 642, 18 
(2006); T. Hambye, K. Kannikc, E. Ma and M. Raidal, Phys. Rev. D 75, 095003 (2007). 

[10] S. L. Glashow and S. Weinberg, Phys. Rev. D 15, 1958 (1977). 

[11] M. Aoki, S. Kancmura, K. Tsumura and K. Yagyu, Phys. Rev. D 80 (2009) 015017. 

[12] V. D. Barger, J. L. Hewett and R. J. N. Phillips, Phys. Rev. D 41, 3421 (1990). 

[13] Y. Grossman, Nucl. Phys. B 426, 355 (1994). 

[14] S. Su and B. Thomas, Phys. Rev. D 79 (2009) 095014; H. E. Logan and D. MacLennan, Phys. Rev. D 
79 (2009) 115022. 

[15] J. F. Gunion, et al., "The Higgs Hunters's Guide" (Addison Wesley, 1990). 

[16] A. Baldini, Nucl. Phys. Proc. Suppl. 168, 334 (2007). 

[17] E. Barberio et al. [Heavy Flavor Averaging Group], arXiv:0808.1297 [hep-cx]. 

[18] G. D. Moore, Phys. Lett. B 439, 357 (1998); Phys. Rev. D 59, 014503 (1998). 

[19] Y. D. Kim, Phys. Atom. Nucl. 69, 1970 (2006); D. S. Akerib, et al., Phys. Rev. Lett. 96, 011302 (2006). 

[20] J. McDonald, Phys. Rev. D 50, 3637 (1994); for a recent study, see e.g., H. Sung Chcon, S. K. Kang 
and C. S. Kim, J. Cosmol. Astropart. Phys. 05 (2008) 004. 

[21] S. Kancmura, Y. Okada and E. Senaha, Phys. Lett. B 606, 361 (2005). 

[22] M. Battaglia, E. Boos and W. M. Yao, |arXiv:hep-ph/0111276| Y. Yasui, et al., |arXiv:hep-ph/0211047| 

[23] S. Kanemura and C. P. Yuan, Phys. Lett. B 530, 188 (2002); Q. H. Cao, S. Kanemura and C. P. Yuan, 
Phys. Rev. D 69, 075008 (2004). 

[24] S. Kancmura, T. Kasai, G. L. Lin, Y. Okada, J. J. Tseng and C. P. Yuan, Phys. Rev. D 64, 053007 
(2001). 

[25] B. K. Bullock, K. Hagiwara and A. D. Martin, Phys. Rev. Lett. 67, 3055 (1991). 



8 



